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ABSTRACT 

The first steps of planet formation are marked by the growth and crystalliza- 
tion of sub-micrometer-sized dust grains accompanied by dust settling toward 
the disk midplane. In this paper we explore whether the first steps of planet 
formation are affected by the presence of medium-separation stellar companions. 
We selected two large samples of disks around single and binary T Tauri stars 
in Taurus that are thought to have only a modest age spread of a few Myr. The 
companions of our binary sample are at projected separations between ~10 and 
450 AU with masses down to about 0.1 Mq. We used the strength and shape of 
the 10 ixm silicate emission feature as a proxy for grain growth and for crystal- 
lization respectively. The degree of dust settling was evaluated from the ratio 
of fluxes at two different mid-infrared wavelengths. We find no statistically sig- 
nificant difference between the distribution of 10 /im silicate emission features 
from single and binary systems. In addition, the distribution of disk flaring is 
indistinguishable between the single and binary system samples. These results 
show that the first steps of planet formation are not affected by the presence of 
a companion at tens of AU. 
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planetary systems: protoplanetary disks 



Introduction 

Two-third of the G stars in the solar 
sighborhood are members of multiple- 



star systems (e.g. IDuquennoy fc Mayor 
199ll ). These binaries and multiple sys- 
tems are often found to harbor gian t 



planets (e.g. iBonavita &: Desideral 120071 ). 



Similarly, young low-mass pre-main se- 
quence stars are very frequently mem- 
bers of multiple systems , mostly bina- 



ries (IMathieu et al.l l2000l : iDuchene et al. 



20071 ). This suggests that planet formation 



around single stars such as our Sun may 
be atypical and urges us to understand the 
effects of stellar companions on planet for- 
mation. We tackle this question from an 
observational point of view. Because nu- 
merical simulations of grain agglomeration 
suggest short timescales for the formation 
of planetesimals (on ly a few 10^ yr, e.g. 
Beckwith et "aDl2000f ). it is crucial to know 



lution is not uniquely controlled by stellar 
age and luminosity but at least one addi- 
tional parameter is present. 

There are two studies suggesting that 
stellar multiplicity could play a major role 



i n the initial dust processing. iMeeus et al. 



how stellar companion (s) affect the dust 
processing in the first few million years. 

Grain growth and the settling of dust 
grains towards the disk midplane are 
thought to represent the first steps in the 



(I2OO3I ) found that among three coeval 
T Tauri disks in the Chamaeleon I star- 
forming region the closest binary sys- 
tem (projected separation of ~120AU) 
sports the strongest contribution from 
large (~ 2 /im) grains and has the high- 
est crystallin e mas s fraction. Similarly, 
Sterzik et all fl2004 h pointed out that the 



disk of a young brown dwarf with a com- 
panion at >30AU shows more processed 
dust than the disk around a single brown 
dwarf. Although the small samples inhibit 
any firm conclusions, these results suggest 



I993I for a review). The study of disks 



around intermediate-mass stars also in- 
dicates a link between grain growth and 
crystallinity. High crystallinity was found 
only when grains larger than the dominant 
sub-micron interstellar g rains were present 
(Ivan Boekel et al.l |2005| ) . In the context 
of these fin dings and dust evolution mod - 



plane t-formation process (e.g. Lissauer &: Stewartl hat companions might trigger rapid dust 

evolution. Intuitively this may happen in 
different ways. A companion could speed 
up dust evolution by dynamically stirring 
the circumstellar dust grains and leading 
to an enh anced collision and gr ain growth 
rate (e.g. iDubrulle et al.lll995l ). In addi- 
tion, the dynamical stirring may lead to 
an increased mixing that could also expose 
larger amounts of dust to temperatures 
high enough (>800K) to be crystallized. 

In this paper we compare two carefully 
constructed samples of disks around single 
and binarjQ stars with a narrow age spread 
to test the hypothesis that binary systems 
have disks with more processed dust and 
flatter structures. In Sect. [2] we describe 
our samples. The data reduction of the 
Spitzer spectra and of the 24 fim MIPS 



els (e.g. iDuUemond fc DominikI 120041 ). 



older disks are expected to have more pro- 
cessed dust (larger grains and crystals) 
than younger disks. In addition, their disk 
structure should be flatter because of the 
gradual settling of large dust grains to- 
wards the disk midplane. However, recent 
observations show that the degree of dust 
processing can be very different even for 
coeval disks around stars of similar spec- 
tral ty pe in the same star-for r aing region 
(e.g. IPrzygodda et al.l l2003l : lApai et al. 



20051 ) . This demonstrates that dust evo- 



few of the systems in our study are triple or 
quadruple systems. For simplicity, we refer to the 
whole sample as binaries) 



2 



photometry is presented in Sect. [31 We 
summarize our results in Sect. H] and dis- 
cuss in Sect. |5] their implications on planet 
formation in single and binary systems. 

2. Sample Definition 

Testing whether stellar companions pro- 
mote the first steps of planet formation 
requires: a) two coeval samples of disks 
around single and binary stars with pre- 
cisely determined multiplicity; b) iden- 
tical spectral type distributions for the 
two samples; c) objects in regions with 
no diffuse PAH emission that could con- 
taminate the dust emission features of the 
targets; d) disks with faint or no PAH 
emission features (for the same reason as 
b). The Taurus-Auriga star-forming re- 
gion has a reasonably complete census of 
its pre-main-sequence stellar population 
and meets best our requirements among 
the nearby star-forming regions. To sat- 
isfy criteria d) we selected disks around 
low-mass stars because they have an or- 
der of magnitude lower PAH emission 
than disks around intermediate-mass stars 



combination of available high-resolution 
imaging and interferometry, spec troscopy, 



(IGeers et al.ll2006l ). The observed age dis- 
tribution of low-mass stars in Taurus can 
be well approximated with a gaussian cen- 
tered at 1.6 Myr and a spread no large r 
than about 2-3 Myr (IHartmannl l200lh . 
This narrow age spread minimizes any pos- 
sible trend of disk evolution with age. In 
summary our sample of disks around sin- 
gle and binary stars is drawn from the 
Taurus-Auriga population of low-mass T 
Tauri stars (TTSs) with circumstellar disks 
and ages between ~l-3 Myr. 

Because it is critically important to in- 
clude only stars with known multiplicity, 
we first selected each target based on a 



1993 


Leinert et al. 


1993 


White & GhezI 


2001). I 



refer to a star as single if it has no known 
companions with brightness above the de- 
tection limit within 10" (or ~1400AU 
at the ~140pc distance of the Taurus- 
Auriga star-forming region). The typical 
detection limit of high-resolution imaging 
surveys is 2-3 mag fainter than the pri- 
mary star in K-band, which corresponds 
to a very-low mass star of ~0.1Mq at 
~2 Myr according to the Baraffe et al. 
1998 isochrones. Only a few binary star 
systems have been observed with tech- 
niques reaching higher contrast (see, e.g. 
the discovery of a brown dwarf co mpan- 
ion to DH Tau bv lltoh et aDl2005l ). The 
typical smallest separation resolvable with 
the imaging surveys is ~0'.'l (~10AU), 
meaning that a few unresolved close bina- 
ries may have been classified as single. Of 
these, binaries with periods < 100 days can 
only modestly contaminate our single star 
sample since their frequenc y in pre-main 



seque nce stars is 8%±3% (IMathieu et al. 



2000l ). Due to the lack of general un- 
derstanding of the separation distribu- 
tion of close pre-main sequence binaries, 
it is not possible to determine the ex- 
act level of contamination. However, the 
required infrared excess (see below) en- 
sures circumstellar dust within a few AU, 
likely inconsistent with the presence of 
a close stellar companion. We also ex- 
cluded a ll known spectroscopic binaries in 
Taurus (IMathieul Il994j : IJensen fc Mathieu 



19971 ) both from the sample of single and 
binary stars. Our binary sample con- 
sists of stars with companions between 
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~0'.'l and 3" thus covering the popula- 
tion of medium-separation stellar com- 
panions at projected separations between 
~14-420AU. After defining the samples, 
we searched for available near- and mid- 



infrared data 


(Kenvon & Hartmann 


1995; 


Stassun et al. 


2001; 


Hartmann et al. 


2005; 


McCabe et al. 


2006; 


Furlan et al. 


20061) to 



select only those targets that have ex- 
cess emission indicative of a circumstel- 
lar disk. These selection criteria resulted 
in a sample of mostly classical TTSs for 
the single and for the primary star in 
the multiple systems. Only CZ Tau and 
IQ Tau in our sar nple are classif i ed as 
weak-lined TT SlI (|Whitefc_Ghd lioOl 



McCabe et al.l 120061 ). The T Tauri types 
for the components of the binary sam- 
ple are mostly classical with the excep- 
tion of FX TauB, IS TauB, V710 TauB, 
and V807 TauB that are weak-lined TTSs 



(Duchene et al. 


1999; 


White & Ghez 


2001; 


Hartigan & Kenvon ' 


2003). Finally, ^ 


Ne se- 



lected only systems with fluxes >0.1Jy 
to ensure high signal-to-noise around the 
10 /im silicate emission feature. 

Our list contains 44 sources, of which 
23 are disks around single stars and 21 
are multiple systems (see Tables [1] and 
[2]). Spectral types for the stellar com- 
panion(s) are available for 18 out of 21 
objects. Most companions have M spec- 
tral type (see Table [2]). The single stars 
(or the primaries of multiple systems) are 
mostly between K5 and M4 spectral type. 
Note also that the distribution of spectral 
types for the single stars is almost iden- 
tical to that of the primary component 
in multiple systems (Fig. [1]). In addition. 



based on their narrow Ha emission lines 



the fairly narrow range in spectral type 
minimizes any difference in the dust com- 
position due to lumi nosity effects. In fact 



the re cent study by iKessler-Silacci et al. 



(120061 ) shows that the shape and strength 
of the 10 //m silicate features of K versus 
M stars are consistent with being drawn 
from the same population. Most of our 
multiple systems are binaries with sepa- 
rations less than 1" (see Table |2]). Sev- 
enteen of them have circumprimary and 
circumsecondary disks based on the K-L, 



try of both components ( 


White & Ghez 


2001 


Hartmann et al. 


|2005; 


McCabe et al. 



20061). DF Tau B and FO Tau B are the 
closest binary components in our sample 
at ~0'.'l. They are also very likely sur- 
rounded by circumstellar disks based on 
accretion signatures like the Ha equiva- 
lent width and the amount of o ptical veil- 
ing (IHartigan fc KenyonI l2003l ). Finally, 



V710 Tau B and V807 Tau B have K- 
N and K-L colors consistent with pho- 
tospheric emission suggestin g that they do 



not have circumstellar disks (White fc Ghez 



20011 ; iMcCabe et al.l 120061 ). A summary 



(jHerbig fc Bel]|ll988D 



of the disk configurations is provided in 
column seven of Table [H Resolved mid- 
infrared photometry is available for 11 ou t 
of 21 binary sources (iMcCabe et al.ll2006l ). 
Another 8 s ources have resolved L-band 
photometry (Iwhite fc Ghea I2OOII ). while 
the two closest binaries FO Tau and DF 
Tau have been resolved only in K-band 
(IWhite fc Ghezll200lh . For at least two- 
third of the binary sample the Spitzer/IKS 
spectrum is dominated by the disk around 
the primary star as indicated by the near- 
and mid-infrared flux ratios of the pri- 
mary and secondary components (column 
eight of Table [2]) . We also note that our 
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Fig. 1. — Histogram of the spectral types for the single and binary systems. The spectral 
type distributions of the single stars and of the primary stars in the binary sample are almost 
identical. 
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selection criteria did not exclude possibly 
more evolved disks such as transition disks. 
GM Aur and DM Tau in our list of single 
stars are amon g the best studied t ransition 
disks (see, e.g. ICalvet et al.l 120051 ). 



3. Data Reduction 

We use the low-resolution Spitzer spec- 
tra around the 10 /im silicate emission fea- 
ture to characterize the degree of dust 
processing in single and binary systems. 
Observations at wavelengths longer than 
~ 20 /im are necessary to trace the evolu- 
tion of the disk structure. For this, we pre- 
fer to use the MIPS 24 /im photometry that 
is available for all (except 6) targets rather 
than the spectra acquired with the long- 
low module of the IRS, which is lacking for 
14 targets. In the following we describe the 
data reduction of the Spitzer IRS spectra 
and of the 24 fim MIPS photometric data. 

3.1. SpitzerlUS low— resolution spec- 
tra 

The IRS data presented in this pa- 
per have been acqui red as part of the 
IRS/GTO program fiFurlan et al.l l2006h 
and became available to the community 
at the end of the year 2005. Only six of 
the objects we selected (CY Tau, DS Tau, 
FM Tau, IS Tau, V710 Tau, and V807 Tau) 
were observed in staring mode, with the 
targets placed in two nod positions along 
the spatial direction of the slit at 1/3 and 
2/3 of the slit length. On all the other 
targets, a 2x3 step mapping observation 
was carried out: the three steps were cho- 
sen to be separated by three-quartes (for 
SL) or half (for LL) of the slit width in 
the dispersion direction and the two steps 
were separated by one-third of the slit in 



the spatial direction. The exposure time 
at each nod position was 6 second for all 
targets with the exception of V710 Tau 
which had a 14 second exposure time per 
nod. 

We downloaded the low-resolution IRS 
data that were processed with the SSC 
pipeline S13.2.0. Our data reduction starts 
from the droopres intermediate data prod- 
uct an d follows the ste p s out lined in de- 
tail in iBouwman et al.l (120071 ). In brief. 



we first subtracted the pairs of imaged 
spectra acquired along the spatial direc- 
tion of the slit in order to correct for 
the background emission and stray-light. 
Then we replaced bad pixels by interpo- 
lating over neighboring, good pixels. Al- 
though for our analysis we use only the 
SL data, we also extracted the LL part of 
the spectra for comparison with the MIPS 
24 /im photometry where both datasets 
are available. Spectra were extracted 
from the background-subtracted pixel- 
corrected images using a 6.0-pixel and 5.0- 
pixel fixed-width aperture in the spatial 
direction for the SL (5.2-14 /im) and LL 
(14 — 35 /im) modules, respectively. The 
low-level fringing at wavelengths >20 /im 
was removed using the irsfringe package 
fiLahuis fc Boogertll2003h . 

Because peak-up images were not ac- 
quired for the majority of the sources, tar- 
gets may not have been positioned accu- 
rately at the center of the slit. We deter- 
mined the actual position of each source 
during extraction by finding the peak emis- 
sion of the wavelength-collapsed source 
profile. Once the spectra are extracted for 
each order, nod and cycle, we computed 
a mean spectrum for each order and as- 
signed as uncertainty at each wavelength 
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Table 1 

Main properties of the selected single TTSs. 



# 


Source 


2MASS J 


Adopted 


Ref 










SpT 


(SpT) 


SI 


AA Tau 


04345542H 


-2428531 


K7 


1, 2 


S2 


BP Tau 


04191583H 


-2906269 


K7 


2, 3 


S3 


CI Tau 


04335200H 


-2250301 


K7 


1, 2 


S4 


CW Tau 


04141700H 


-2810578 


K3 


1, 2 


S5 


CX Tau 


04144786H 


-2648110 


M2.5 


1, 2 


S6 


CY Tau 


04173372H 


-2820468 


M2 


2, 3 


S7 


DE Tau 


04215563H 


-2755060 


Ml 


2, 3 


S8 


DH Tau* 


04294155H 


-2632582 


M2 


2, 4 


S9 


DL Tau 


04333906+2520382 


K7 


1, 2 


SIO 


DM Tau* 


0433487H 


-1810099 


Ml 


1, 2 


Sll 


DN Tau 


04352737H 


-2414589 


MO 


1, 2 


S12 


DO Tau 


04382858H 


-2610494 


MO 


1, 2 


S13 


DP Tau 


04423769H 


-2515374 


M0.5 


1, 2 


S14 


DS Tau 


04474859H 


-2925112 


K5 


1, 2 


S15 


FM Tau 


04141358H 


-2812492 


Ml 


2, 4 


S16 


FN Tau 


04141458H 


-2827580 


M5 


2 


S17 


FZ Tau* 


04323176H 


-2420029 


MO 


2, 4, 5 


S18 


Gl Tau* 


04333405H 


-2421170 


K7 


2, 4 


S19 


GK Tau* 


04333456H 


-2421058 


K7 


2, 4 


S20 


GM Aur* 


04551098H 


-3021595 


K3 


1, 2 


S21 


IP Tau 


04245708H 


-2711565 


MO 


1, 2 


S22 


IQ Tau 


04295156H 


-2606448 


M0.5 


1, 2 


S23 


LkCa 15 


04391779H 


-2221034 


K5 


1, 2 



*These st ars have companions at projected separations just 
above 10" ijHartigan et al.l Il994l ) : DH Tau/DI Tau (15"); 
FY Tau/FZ Tau (16'.'9); GK Tau/GI Tau (12'.'9) 

* Known transition disks (see, e.g. ICalvet et aLllSoOSh . 

Note. — The 2MASS source name includes the J2000 sexa- 
gesimal, equatoria l position in the form: hhmmssss-|-ddmmsss 
ijCutri et al.ll2003h . 

References. ~~ (1) iHerbig fc Belil 11988"; (2) 

kcnvon fc HartmannI Il995l: (3')[strom fc Strom ,1994 : (4) 
Hartigan et al.lll994l: (5') IWhite fc Ghej|200lf 



7 



Table 2 

Main properties of the selected binary TTSs. 



# 


Source 


2MASS J 


Adopted 


Separation 


Ref 


Disk 


Flux ratio 


Ref 










SpT* 


(arcscc) 


(SpT) 


configuration* 


(Alter)* 


(flux ratio) 


Bl 


CoKu Tau3 A-B 


04354093H 


-2411087 


Ml 


2.05 


1 


cp+cs 


4.7 (L) 


3 


B2 


CZ Tau A-B 


04183158H 


-2816585 


Ml. 5 


0.32 


1 


cp+cs 


1.71 (L) 


10 


B3 


DD Tau A-B 


04183112H 


-2816290 


M3+M3 


0.56 


2, 3 


cp+cs 


1.75 (N) 


10 


B4 


DF Tau A -B 


04270280H 


-2542223 


M0.5+M3 


0.09 


4, 3 


cp+cs 


1.6 (K) 


3 


B5 


DK Tau A-B 


04304425H 


-2601244 


K9+M1 


2.30 


1, 5 


cp+cs 


8.53 (SiC) 


10 


B6 


FO Tau A-B 


04144928H 


-2812305 


M2+M2 


0.15 


4, 3 


cp+cs 


1.7 (L) 


3 


B7 


FS Tau A-B 


04220217H 


-2657304 


M1+M4 


0.23 


4, 3 


cp+cs 


5 (L) 


3 


B8 


FV Tau A-B 


04265352H 


-2606543 


K5+K6 


0.72 


4, 2, 3 


cp+cs 


2.2 (N) 


10 


B9 


FX Tau A-B 


04302961+2426450 


M1+M4 


0.89 


1, 5 


cp+cs 


2.8 (N) 


10 


BIO 


GG Tau Aa-Ab 


04323034+1731406 


K7+M0.5 


0.25 


6 


cp+cs 


1.03 (N) 


10 


Bll 


GH Tau A-B 


04330622H 


-2409339 


M1.5+M2 


0.31 


3 


cp+cs 


1.45 (N) 


10 


B12 


GN Tau A-B 


04392090H 


-2545021 


M2.5 


0.33 


7 


cp+cs 


1.53 (L) 


10 


BIS 


HK Tau A-B 


04315056H 


-2424180 


M1+M2 


2.34 


1, 5 


cp+cs 


30 (SiC) 


10 


B14 


HN Tau A-B 


04333935H 


-1751523 


K5+M4 


3.11 


2 


cp+cs 


65 (L) 


3 


B15 


IS Tau A-B 


04333678H 


-2609492 


K7+M4.5 


0.22 


8, 3 


cp+cs 


9 (L) 


3 


B16 


IT Tau A-B 


04335470H 


-2613275 


K3+M4 


2.39 


5 


cp+cs 


2.95 (SiC) 


10 


B17 


RW Aur A-B 


05074953H 


-3024050 


K1+K5 


1.42 


9, 5 


cp+cs 


13.63 (N) 


10 


B18 


UY Aur A-B 


04514737H 


-3047134 


K7+M2 


0.88 


1, 5 


cp+cs 


2.07 (N) 


10 


B19 


V710 Tau A-B 


04315779H 


-1821380 


M0.5+M2 


3.17 


2 


cp 


12.7 (SiC) 


10 


B20 


V807 Tau A-B 


04330664H 


-2409549 


K7+M3 


0.3 


4, 2, 3 


cp 


3.6 (L) 


3 


B21 


V955 Tau A-B 


04420777H 


-2523118 


K5+M1 


0.33 


3 


cp+cs 


5.6 (L) 


3 



*The second spectral type, when available, is for the secondary star. 

*'cp' stands for circumprimary disk while 'cs' stands for circumsecondary disk. These disk conflgurations have been determined 
via resolved optical and infrared photometry, see Sect.[2]for details. 

*Flux ratios are calculated as prim ary /secondary. In pa renthesis we provide the fllter at which the flux ratio has been calculated 
(the SiC filter is centered at 11.8 um. iMcCabe et aI]|2006D . 



Note. — The 2 MASS source name includes the J2000 sexagesimal, equatorial position in the form: hhmmssss+ddmmsss 
l lCutri et al.ll2003h . 

Refere nces. — (DiLeinert et al.ll993: ('2 VHartigan et al.'l994': (3) 'White & Ghcz"200l'; ('4VCohe n_fc Kuhill979l: (SllDuchene et al, 

Il999l : (e'l lWhite et al.lll999l : (7') IWhite fc Basri 2003 : (8) Martin ct al. 1994 : (9) Mundt &: Giampa'palll982l : (lOl lMcCabe et al.ll200e 
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the 1-sigma standard deviation of the dis- 
tribution of the data points. 

The absolute flux calibration was done 
using order-based spectral response func- 
tions created within the Formation and 



Spitzer Legacy program ( 


Hines et al 


2005 


Mever et al. 


2006; 


Bouwman et al. 


2007) 



The advantages of this spectral response 
function over the standard SSC bed cali- 
bration are: a) the use of a larger number 
of calibrators from the FEPS program; b) 
the spectral response function is order and 
nod-position based; c) it is a ID spectral 
response function allowing a better rejec- 
tion of bad pixels. The estimated abso- 
lute fl ux calibration uncerta inty is around 



10% (iBouwman et al.l 120071 ) and is prop- 



agated to the flux uncertainties assigned 
at each wavelength. The resulting cali- 
brated spectra are available in the elec- 
tronic edition of the Astrophysical Journal 
(Figs. El m [ini [m [m [131 [H and [ISD . 

We checked the absolute flux calibration 
of the IRS SL module usi ng the published 
IRAC 8 fim photometry by 



( 120051 ) and iLuhman et al. 



Hartmann et al. 



(120061 ). whose 



data have been acquired in two differ- 
ent campaigns, February-March 2004 and 
February 2005 respectively. Eleven sources 
have IRAC 8 /im magnitudes from both 
campaigns, while 5 more sources have 
only magnitudes from Fe bruary-March 
2004 (IHartmann et al.ll2005l ). an d 14 more 



sourc es only from February 2005 (iLuhman et al. 
20061 ). The mean difference in the IRAC 



8 nm fluxes for the 1 1 sources observed in 
both campaigns is 16%, much larger than 
the IRAC absolute flux calibration accu- 



racy of a few percent (IReach et al.ll2005l ). 
The largest deviations are for FV Tau, 



DL Tau, and DO Tau whose 8 fim IRAC 
fluxes decreased by factors of 33%, 23%, 
and 23% respectively in a year baseline. 
These differences are likely due to intrin- 
sic stellar variability. The IRS 8 /im fluxes 
integrated over the IRAC 8 fim spectral 
response curve agree on average within 
~10% of the published IRAC photometry, 
which is within the estimated IRS flux cal- 
ibration uncertainty. We also verifled that 
our spect ra very well agree w ith those pub- 
lished by iFurlan et al.l (120061 ) who adopted 
a different data reduction. 

3.2. Spitzer 24 /xm photometry 

MIPS 24 /im (hereafter MIPS24) data 
are available in the Spitzer archive for all 
but 6 of the targets we selected (specifl- 
cally DM Tau, DS Tau, GG Tau, HN Tau, 
RW Aur, and UY Aur). The majority 
of the MIPS data have been acquired as 
part of the Taurus Spitzer Legacy Program 
(PI D. Padgett) in February-March 2005 
(id=3584) and March 2007 (id=30816) us- 
ing the MIPS scan map operational mode. 
V710 Tau and GM Aur are not covered 
by these programs but have data from the 
c2d Spitzer Legacy Program (PI N. Evans, 
pid=173) and the GTO program respec- 
tively (PI G. Fazio, pid=37). We down- 
loaded the post-bed MIPS 24 fim prod- 
ucts processed through the SSC pipelines 
S14.4.0 or later. In the case of photometry 
mode (only for V710 Tau), the SSC prod- 
uct is an averaged and registered single im- 
age while in the case of scan maps (for all 
other sources), the product is a distortion- 
corrected mosaic imag^. The 24 /im post- 
bed images generated from pipeline ver- 



^http://ssc. spitzer. caltech.edu/mips/dh/niipsdatahandbook3. 2.pdf 
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Table 3 

Fluxes at 24 /im from MIPS images (third column, F24) or at 25 ^m from 

IRS SPECTRA (fourth COLUMN, F25). THE ABSOLUTE PHOTOMETRIC UNCERTAINTY 
IS EXPECTED TO BE ~10%. THIS UNCERTAINTY INCLUDES BOTH THE INTERNAL 
RANDOM AND THE ABSOLUTE CALIBRATION UNCERTAINTY, SEE TEXT FOR DETAILS. 



# 


Source 


F24 


F25 


Comment 






[mJy] 


[mJy] 





SI 


AA Tau 


525 






S2 


BP Tau 


666 






S3 


CI Tau 


997 






S4 


CW Tau 


1445 






S5 


C;X Tau 


296 






S6 


CY Tau 


123 






S7 


DE Tau 


677 






S8 


DH Tau 


316 






S9 


DL Tau 


965 






SIO 


DM Tau 




320 


no MIPS24 


Sll 


DN Tau 


424 






S12 


DO Tau 




2951 


MIPS24 saturated 


S13 


DP Tau 


1281 






S14 


DS Tau 




303 


no MIPS24 


SIB 


FM Tau 


462 






S16 


FN Tau 


1047 






S17 


FZ Tau 


1057 






818 


GI Tau 


1006 






S19 


GK Tau 




1699 


MIPS24 saturated 


S20 


GM Aur 


746 






S21 


IP Tau 


277 






S22 


IQ Tau 


544 






S23 


LkCa 15 


398 






Bl 


CoKu Tau3 


327 






B2 


CZ Tau 


1232 






B3 


DD Tau 


1575 






B4 


DF Tau 


945 






B5 


DK Tau 


1286 






B6 


FO Tau 


529 






B7 


FS Tau 




2112 


MIPS24 saturated 


B8 


FV Tau 




2344 


MIPS24 saturated 


B9 


FX Tau 


420 






BIO 


GG Tau 




1266 


no MIPS24 


Bll 


GH Tau 


379 






B12 


GN Tau 


536 






B13 


HK Tau 


823 






B14 


HN Tau 




2988 


no MIPS24 


B15 


IS Tau 


228 






B16 


IT Tau 


264 






B17 


RW Aur 




2012 


no MIPS24 


BIB 


UY Aur 




5918 


no MIPS24 


B19 


V710 Tau 


245 






B20 


V807 Tau 


476 






B21 


V956 Tau 


546 
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sions later than S14 have improved flat 
field corrections and are suitable for pho- 
tometry at an accuracy of 10% (Sect. 9.1 
of the MIPS Data Handbook v 3.2, link in 
footnote 3), sufficient for the purposes of 
this study. 

Ap erture photometry was d one using 
IDP3 fISchneider fc Stobiell2002h . The cen- 
teroid of the aperture was found by fit- 
ting a Gaussian to each source. We used 
an aperture radius of 6'.'6 (2.7 pixels) and 
background annulus from 20" to 32". We 
opted for this intermediate aperture radius 
among those suggested by the Spitzer Sci- 
ence Center to include the emission from 
both binaries (as in the IRS spectra) and 
to exclude the emission from companions 
at > 10" from some of the single stars (see, 
note to Table [1]). To recover the total 
flux we then applied the proper aperture 
correction of 1.648 0. DO Tau, GK Tau, 
FS Tau, and FV Tau have 24 /im fluxes 
>1.6Jy and are thus saturated in the 
3 sec exposures acquired within the Tau- 
rus Spitzer Legacy Program. For these 
sources as well as for the 6 targets that 
lack MIPS24 data, we computed the 25 /xm 
flux (integrated flux from 23.5 to 26.5 /xm) 
from our SL 6 fim flux and t he nfi,^^ spec- 
tral in dices from Table 4 of iFurlan et al. 
fl2006f lFI. We will show in Sect. lU that 



our results are not affected by the use of 
25 fim fluxes instead of 24 /im fluxes for 
these 10 objects. For the sources that 
have both MIPS24 observations and LL 
data we flnd that our aperture photom- 



''http://ssc.spitzer.caltech.edu/mips/apercorr/ 
^most of these sources are so bright that were ob- 
served only with the IRS high-resolution module. 
Therefore we could not compute 24 /im fluxes from 
the LL data we reduced 



etry agrees within about 10% with the 
IRS 24 iim flux integrated over the MIPS24 
spectral response curve. Table [3] summa- 
rize the MIPS24 or the IRS25 fluxes for 
our sample of single and binary systems. 

4. Results 

4.1. Strength and Shape of the 
10 /im Feature 

The 10 yum silicate emission feature 
traces silicate dust in the optically thin 
layer of circumstellar disks out to about 
1 AU from a star with solar luminosity (e.g. 



Kessler-Silacci et al.l 120071 ) . Because the 



smaller the grain the longer its timescale 
to settle to the disk midplane, the 10 fim 
feature probes the population of the small- 
est grains in a disk, up to a few micron in 
size. The strength and shape of the fea- 
ture bear information on the amount of 
processing that dust grains have under- 
gone in the disk. In particular, the flux 
ratio of normalized spectra at 11.3 (or 8.6) 
over 9.8 fim can be used as a proxy for the 
degree of crystallization while the peak-to- 



Bouwman et al 



grain growth (e.g . 

Boekel et a h '20051: lApai et al 



van 





2001 




2005 



Bouwman et al..,2007 ). We use these band 
strengths to search for differences in the 
dust processing of single and binary sys- 
tems. 

We flrst flt a third-order polynomial^ 
to the spectral data outside the 10 /xm sili- 



third order polynomial gives smaller residuals in 
the spectrum minus fitted continuum than a first, 
a second, or a fourth order polynomial. We also 
verified that the results discussed in this section 
do not change when using lower or higher order 
polynomials to fit the continuum. 
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Fig. 2. — Strengths and shapes of the 10 /xm sihcate emission features for the sample of 
single (red circles) and binary (blue squares) systems. On the y-axis we plot the ratios 
of the normalized flux at 8.6 and 11.3 /im over 9.8 /im, which are proxies for the degree of 
crystallinity. The x-axis gives the peak over continuum in the 10 /im region of normalized 
spectra, which indicates the amount of grain growth. Objects with processed dust are located 
on the top left of the diagram. Average error bars are shown on the right upper corner of 
each plot (see Sect. 14. II for details). 
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cate emission feature (between 6 and 8 /xm 
and between 12 and 14 //m) and then nor- 
malize our spectra to the fitted continuum. 
This normahzation ensures that the shape 
of the spectral features remains identical 



to th e original one (e.g. Ivan Boekel et al. 
20051 ) . Fig. [2] shows the band strengths for 
our sample of single (red circles) and mul- 
tiple (blue squares) systems. The flux den- 
sities at 8.6, 9.8, and 11.3 fim are the mean 
values of flux densities within ±0 . 1 fim of 
these wavelengths. To analyze the un- 
certainties we followed a Monte Carlo ap- 
proach. For each disk we added a normally 
distributed noise to the spectrum and com- 
puted the peak-over-continuum and flux 
ratios 500 times. Their standard devia- 
tion gives the uncertainty on the peak- 
over-continuum and flux ratios for each 
target. The error bars in Fig. [2] are the av- 
erage uncertainties. The simulated noise 
had two components: i) a normally dis- 
tributed noise at each wavelength with an 
amplitude equal to the flux uncertainty at 
that wavelength; ii) a random calibration 
uncertainty (equal at all wavelengths) with 
an amplitude of 10%. 

Objects on the bottom right side of 
Fig. [2] have 10 /im emission features sim- 
ilar to the ISM absorption feature, thus 
their dust grains have experienced little 
processing. Conversely, objects on the top 
left side of the plot have disks with larger 
grains (micron in size) and crystals. There 
is an obvious anticorrelation between the 
strength of the 10 fim silicate emission fea- 
ture and the F11.3/F9.8 (or F8.6/F9.8) 
flux ratios, which is apparent not only in 
Fig. [2] but also in other similar plots in 



the literature (e .g. Ivan Boekel et al.l 12005 



Apai et al.ll2005l ). We use the Kendall's r 



test to measure the degree of anticorrela- 
tion. This nonparametric test uses the rel- 
ative rank ordering of p airs of data to corn - 
pute two values (see e.g. IPress et al.lll993l ): 
the rank correlation coefficient r which 
runs between -1 (complete rank reversal) 
and 1 (complete rank agreement), and the 
two-sided probability P that the variables 
are uncorrelated (r = 0). For the single 
stars in Fig. [2] (circles), Kendall's r = —0.8 
and P = 6 X 10"^ for the F8.6/F9.8 flux ra- 
tio versus the peak-over-continuum while 
r = -0.6 and P = 2 X 10"^ for the 
F11.3/F9.8 flux ratio versus the peak- 
over-continuum. For the binary stars in 
Fig. [2] (squares), Kendall's r = —0.5 and 
P < 0.001 for the F8.6/F9.8 and the 
F11.3/F9.8 flux ratios versus the peak- 
over-continuum. We interpret these trends 
as conflrmation that there is a signifl- 
cant anticorrelation between the strength 
of the silicate emission features and the 
flux ratios presented in Fig. [2l There- 
fore, both the peak-over-continuum and 
the F11.3/F9.8 (or F8.6/F9.8) flux ratios 
provide an analogue measurement of dust 
processing in circumstellar disks. 

The main question we want to answer 
in this paper is whether the population of 
disks around single stars statistically differ 
in terms of dust processing from the pop- 
ulation of disks around binaries. In other 
words, do either of the two populations 
have more processed dust in their disks 
than the other? To answer this question we 
apply the Kolmogor ov-Smirnoy (here after 
K-S) test (see, e.g. iPress et al.lll993l ) and 
the Mann-Wh itney U (hereafter MWU) 
test (see, e.g. iRonald fc Raymond! Il985l ) 
on the peak-over-continuum and the flux 
ratios (F11.3/F9.8 and F8.6/F9.8) pre- 
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Fig. 3. — Histogram of the 10 fim feature strengths for the single (red) and binary (blue) 
systems. The K-S test indicates that the two distributions are consistent with having been 
drawn from the same parent population. This suggests that the growth of grains up to a few 
micron in size is not affected by the presence of a medium-separation stellar companion. 
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Fig. 4. — Histogram of the flux ratios at 11.3 over 9.8 /im for the single (red) and binary (blue) 
systems. The K-S test does not indicate that the two distributions differ statistically. This 
suggests that crystalline processing is not influenced by the presence of stellar companions 
at tens of AU. 
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Fig. 5. — Strength (squares) and shape (stars) of the 10 /im sihcate emission feature versus 
projected binary separation. The peak-over-continuum compared to the projected separa- 
tion (squares) gives a Kendall's r=0.1 and P=0.4 suggesting that these variables are not 
correlated. Similarly, the F11.3/F9.8 flux ratio is uncorrelated with the binary projected 
separation (stars) as indicated by a Kendall's r=0.02 and P=0.9. 
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sented in Fig. H The K-S and the MWU 
tests are both non-parametric tests with 
the null hypothesis that two samples are 
drawn from a single population and work 
on unbinned data. Because the K-S test 
gives large probabilities of 0.2 for the peak- 
over-continuum of single and binaries and 
0.3 for their F11.3/F9.8 (and F8.6/F9.8) 
flux ratios, we conclude that both samples 
could have been drawn from the same dis- 
tribution. Similarly the MWU test does 
not reveal any statistically significant dif- 
ference in the degree of grain growth and 
crystallinity of the single and binary sam- 
ples: the probabihty that both distribu- 
tions were drawn from the same parent 
population for the peak-over continuum is 
0.08 while for the flux ratios F11.3/F9.8 
(and F8.6/F9.8) is 0.1. Figs. [3] and ll pro- 
vide a visualization of the peak-over con- 
tinuum and F11.3/F9.8 flux ratio distribu- 
tions for the samples of disks around single 
and binary stars. Further supporting our 
conclusion, Fig. [S] shows no clear trend of 
the strength (or shape) of the 10 fim fea- 
ture with the projected separation of the 
companion. Kendall's r are close to and 
P are large for the rank ordering of peak- 
over-continuum (or F11.3/F9.8) compared 
to the projected separation (see caption of 
Fig. [5]) suggesting that the variables are 
uncorrelated. 

4.2. Dust settling 

Coagulation models show that if the 
dust grains in the upper layers of a flared 
disk become sufficiently large, they will 
gravitationally settle towards the mid- 



a way to search for grain growth is to 
evaluate the flaring of circumstellar disks. 
While the shape of the 10 fim feature is 
sensitive to the presence of grains of a few 
microns in size in the disk upper layer 
the disk flaring should probe the over- 
all grain population of large r grains (e.g. 
DuUemond &: DominikI I2004J ). To evalu- 



disk geometry (e.g. 


SchraDler & Hennine 


2004: 


Nomura & Nakagawa 


,20061). Thus 



ate the disk flaring we use the ratio of 
fluxes at two different mid-infrared wave- 
lengths. This procedure is justifled by two 
facts: a) the continuum flux emitted from 
the surface layer of the disk becomes pro- 
portional to the disk flaring at radial dis- 
tances from solar-type stars of .4 AU or 
larger (IChiang fc GoldreichI 119971 ) and b) 
these radial distances are probed by mid- 
infrared observations. 

We integrate the flux of our SL infrared 
spectra in two wavelength bands one short- 
ward and one longward of the 10 fxm sili- 
cate emission feature: 5.4-6.0 fim (central 
wavelength 5.7/im) and 12.5-14.0 /im (cen- 
tral wavelength 13.25 /im). In addition, 
we use the MIPS24 photometry (or the 
IRS25 flux when MIPS is not available, 
see Table [3]) to trace the disk flaring out 
to a few AU from the central star. The 
calculated flux ratios are shown in Fig. [6] 
and histograms are presented in Fig. [7l 
Larger flaring is indicated by higher ratios 
of long-wavelength continuum to short- 
wavelength continuum flux from the dust 
disk. These plots show that both single 
and binary systems have a large variety 
of disk structures with no preference for 
any structure in the two samples. The 
K-S and the MWU probabilities that the 
the flux ratio distributions for single and 
binary stars come from the same parent 
population are higher than 0.1 conflrming 
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3.0 



Fig. 6. — Flux ratios at 13.25 over 5.7/im (top), at 24.0 over 13.25 //m (middle), and at 24.0 
over 5.7/im (bottom) versus the peak-over-continuum in the 10 /im feature. These ratios are 
a proxy for the disk flaring; more flaring is indicated by larger ratios in the figure. Squares 
represent binaries, while circles indicate single stars. Filled symbols are for those sources 
that do not have MIPS24 photometry or are saturated in the MIPS exposures. For these 
sources we plot the IRS25 flux computed as described in Sect. 13. 2[ Both single and binary 
systems have a large variety of disk structures. 
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Fig. 7. — Histograms for the flux ratios at 13.25 over 5.7/im (top), at 24.0 over 13.25 /im 
(middle), and at 24.0 over 5.7//m (bottom) for the single (red) and binary (blue) systems. 
The distributions of flux ratios (disk flaring) do not differ statistically according to the K-S 
and the MWU tests. MWU probabilities are 0.22, 0.13, and 0.18 for the F13/F6, F24/F13, 
and F24/F6 respectively. 
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that the distributions of disk flaring for the 
two samples are not statistically different 
(the same result holds when excluding the 
sources with IRS25 fluxes). 

5. Discussion 

Our study shows a large diversity in the 
10 /im silicate emission features and SED 
slopes of T Tauri disks. We found that 
neither the dust processing nor the disk 
flaring correlates with the multiplicity of 
the sources. These results are particularly 
interesting for two aspects that will be dis- 
cussed in the following. 

5.1. Medium— separation binaries and 
planet formation 

A stellar companion induces tidal forces 
in a disk that become particularly strong 
at resonance points. Resonant interac- 
tions result in the excitation of density 
waves that can truncate a disk and act 
to modify the binary eccent ricity (see e.g. 
Lubow fc Artymowicz 2000l for a review). 
Theoretical calculations of binary-disk in- 
teractions predict that circumstellar disks 
will be truncated at 0.2-0.5 times the bi- 
nary semimajor axis a, with the exact val- 
ues depending on e ccentricity, mass ratio, 
and disk viscosity (lArtymowicz fc Lubow 
I994J ). These theoretical expectations are 
supported by millimeter observations of bi- 
naries tracing the optically thin dust emis- 
sion and thus the total disk mass in the 
system. There is evidence for a dimin- 
ished millimeter flux (hence disk mass) 
among the 1-100 AU binaries in compar- 
ison to wider binari es or single stars (e.g. 
Mathieu et al.ll2000l for a review). This re- 
sult is qualitatively consistent with the cir- 
cumstellar disks of medium-separation bi- 



naries being tidally truncated at 0.2-0. 5a. 
Two-thirds of our sources have stellar com- 
panions between 0.1"-1", with a mean pro- 
jected separation of 0.4" or 56 AU at the 
distance of Taurus. Therefore, the typical 
truncation radius for disks in our sample 
is >11-28AU, well outside the location of 
Jupiter and Saturn in our Solar System. 
Even in other systems these outer radii 
are f ound to be devoid of giant planets 
(e.g. iKasper et al.l 120071 ). This fact sug- 



gests that the formation of terrestrial and 
giant planets may proceed undisturbed in 
disks around medium-separation binaries 
even if these disks are constrained in size. 

Early investigations of young TTSs 
found no signiflcant difference in the fre- 
quency of near- and mid-infrared ex- 
cess emission betwe en single and binary 



star systems (e.g. ISimon fc Pratd Il995 



Jensen et al.lll996l l With the 60^m IRAS 



flux probing dust <10AU from the cen- 
tral star, these measurements demonstrate 
that binary systems as ofte n have disks 



as sin gle stars do. Recently iMonin et al. 



(I2OO7I ) analyzed the separation distribu- 
tions of binaries with and without disks 
and found no statistical difference. Since 
most of their binaries have projected sep- 
arations >20AU, their result shows that 
medium-and wide- separation binaries do 
not have a signiflcant effect on the cir- 
cumstellar disk lifetime. Our work in- 
dicates that these disks also evolve in a 
similar way. The extent of dust process- 
ing in the disk surface layer and the de- 
gree of dust settling in binary disk sys- 
tems do not statistically differ from those 
in disks around single stars. This sug- 
gests that the flrst few Myr of disk evo- 
lution in the terrestrial (and maybe out 
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to the giant) planet-forming region are 
not affected by medium-separation stel- 
lar companions. Whether the disk evo- 
lution proceeds undisturbed for tens of 
millions of years until planets are fully 
formed c annot yet be a s sessed observa- 
tionally. iBouwman et al.l (120061 ) estimate 
a mean disk dispersion timescale of ~5 Myr 
for close (<4AU) binaries in contrast to a 
timescale of ~9 Myr for single star systems. 
They argue that the time available to form 
planets in close binary systems is consid- 
erably shorter than that in disks around 
single stars, which may inhibit planet for- 
mation. The only two medium-separation 
binaries in their sample hint for a disk 
dispersal timescale comparable to that of 
single stars suggesting a similar disk evo- 
lution for single and medium-separation 
binary systems over the first ~ 10 Myr. 

Exoplanet surveys offer us a glimpse 
into the frequency and properties of giant 
planets in multiple st a r syst ems. Recently 
Eggenberger fc Udryl (120071 ) reported 42 
planets orbiting bina ry and multiple stars 
(see, t heir Table 1). iBonavita fc Desidera 
(120071 ) analyze a subsample of radial veloc- 
ity planet host stars with uniform planet 
detectability and demonstrate that the 
overall frequency of giant planets in bi- 
naries is not statistically different from 
that of planets in single stars. However, 
they find indications for a lower frequency 
of radial velocity planets in the subgroup 
of close- and medium-separation bina- 
ries (< 50 — 100 AU). In a corn pl emen- 



tary study, iDesidera fc Barbieril (120071 ) 



find that the mass distribution of planets in 
binaries with separations < 300 — 500 AU 
is statistically different from that around 
wider binaries and single stars: Massive 



planets in short-period orbits are found 
predominantly around close- and medium- 
separation binaries. Taken together, the 
results from the frequency and properties 
of exoplanets suggest that a stellar com- 
panion with separation less than a few 
hundred AU affects giant planet formation 
and/or the subsequent migration. Numer- 
ical si mulat i ons s eem to support this no- 
tion. iKley (l2000h shows that a fairly ec- 
centric (cbin = 0.5) stellar companion at 
50-100 AU enhances the growth rate of a 
Jupiter mass planet embedded in a circum- 
stellar disk and makes its inward migra- 
tion m ore rapid. Recently, iKley fc Nelson 
(120071 ) confirm these trends by following 
the evolution of a 30 M^ protoplanet in a 
disk truncated by a stellar companion at 
18.5 AU and Cbm = 0.36, like the 7 Cep 
binary system. Our study shows that the 
early evolution of protoplanetary disks sur- 
rounding binary stars is similar to that in 
single stars indicating that that the differ- 
ences in the exoplanet properties arise in 
the later stages of their formation and/or 
migration. 

Whether terrestrial planet formation is 
also affected by medium-separation bina- 
ries cannot be yet addressed observation- 
ally. Our study shows that the initial 
dust processing is not impacted by the 
presence of a stellar companion. Based 
on the fact that the build-up of plan- 
etesimals as large as the ~500-km Vesta 
has occurred in th e first 3.8=bl . 3 My r of 
the Solar nebula (IKleine et al.l |2002| ). it 
is reasonable to speculate on the basis of 
our study that the formation of planetesi- 
mals in binary and single systems proceed 
along, if not on identical avenues. An- 
other indication supporting this sugges- 
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tion comes from the finding of a similar 
incidence of debris d isks in Gyr-old sin - 
gle and binary stars (ITrilling et al.l 120061 ). 
If the debris dust is produced by colhd- 
ing asteroids, then the similar rate of de- 
bris dust in binaries implies that plan- 
etesimal formation is not inhibited by the 
presence of stellar companions. Recent 
simulations of the later stages of terres- 
trial planet formation show that rocky 
planets can f orm in a wide variety of b i- 
nary systems (iQuintana &: Lissauerl 120071 ). 
The binary periastron is the most impor- 
tant parameter in limiting the number of 
forming planets and t heir r ange of orbits. 



Quintana fc Lissauerl (120071 ) show that bi 



naries with periastron >10 AU, comprising 
most of the medium-separation binaries 
investigated in this paper, can form ter- 
restrial planets over the entire range of 
orbits allowed for single stars. As a re- 
sult more than 50% of the binary systems 
in th e Milky Way (jPuquennoy fc Mayor 



199ll ) are wide enough to allow the forma- 



tion of Earth-like planets. 

5.2. The diversity in silicate fea- 
tures and SEDs 

Although small sample statistics sug- 
gested a correlation between stell ar multi- 
plicity and initial dust processing fjMeeus et al, 



2003; Sterzik et al.ll2004J : ISiciha-Aguilar et al. 



20071), our study demonstrates that medium- 



separation stellar companions do not ap- 
preciably affect the growth and crystalliza- 
tion of dust in circumstellar disks. Given 
the criteria applied to select our samples, 
we can also exclude that age, spectral type, 
and stellar environment can account for 
the large variety of observed silicate emis- 
sion features and SED slopes in our study. 



There may be several other factors con- 
tributing to this diversity that will be fully 
explored in an upcoming contribution. In 
the following we briefly mention two of 
them: 

Turbulence in circumstellar disks not 
only drives the accretion of gas onto the 
central star but also replenishes the disk 
atmosphere with more grains that can be 
larger in size. If the grains inferred from 
the 10 /im silicate emission feature reflect 
the level of disk turbulence, the strength 
of the features shoul d depend on the stel 



l ar ac cretion rates. ISicilia-Aguilar et al. 



(120071 ) note that stars with strong features 
tend to have large accretion rates in their 
sample of several Myr old intermediate- 
and low-mass stars. This trend may be 
the result of turbulence determining the 
grain population in the disk atmosphere. 
Alternatively, the trend could be due to 
the more massive stars (that have typi- 
cally larger accretion rates) in their sam- 
ple heating larger disk area and thus pro- 
ducing st ronger silicate emission fe atures 
(see, e.g. iKessler-Silacci et al.l 120071 ). The 
ten tative correlation se e n in the sample 
of ISicilia-Aguilar et al.l (120071 ) needs to 



be confirmed using a larger and more 
homogeneous sample of stars with well- 
determined accretion rates. 

Different initial conditions for the col- 
lapsing cores may also leave their imprints 
on the formation and evolution of circum- 
stellar dis ks. This possibility has b een ex- 
plored by iDuUemond et al.l (120061 ) to ex- 
plain crystallization of dust grains in the 
early stages of disk evolution. In their 
model the level of crystallinity depends 
crucially on the rotation rate of the collaps- 
ing cloud core because this determines the 
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radius at which the infaUing matter reaches 
the disk: rapidly rotating clouds would 
evolve into disks with low crystallinity, 
while slowly rotating clouds into disks with 
high crystallinity. 

6. Summary 

In this paper we explored the effect of 
a stellar companion on the initial growth 
and settling of dust grains in circumstellar 
disks. We constructed two large samples 
of disks around single and binary TTSs 
with a narrow age spread and a spectral 
type distribution for the single stars iden- 
tical to that of the primary stars in the 
binary sample. We used the strength of 
the 10 /xm silicate emission feature derived 
from IRS/ Spitzer spectra as a proxy for 
grain growth and the SED slope of circum- 
stellar disks as a proxy for dust settling. 
Our results can be summarized as follows: 

- there is no statistically significant differ- 
ence between the distribution of 10 /xm sil- 
icate emission features from single and bi- 
nary systems. 

- the distribution of disk flaring is indis- 
tinguishable between the single and binary 
system samples. 

These results show that stellar compan- 
ions at projected separations of > 10 AU do 
not appreciably affect the degree of crys- 
tallinity nor the degree of grain growth. 
Based on the combination of these and 
other results we argue that the formation 
of planetesimals and possibly terrestrial 
planets is not inhibited in a circumstel- 
lar disk perturbed by a medium-separation 
stellar companion. 

This work is based on observations 
made with the Spitzer Space Telescope, 
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